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Abstract. We compute the shape derivative of the functional Q — » Ai q, 
where Ai q denotes the first eigenvalue of the 1-Laplacian on fl. As an appli- 
cation, we find that the ball is critical among the volume-preserving deforma- 
tions. 



Let 51 be a smooth bounded domain in K™, n > 2. The 1-Laplacian on SI is the 
formal operator 



Aiu = — div 



|Vu| 

we get by a formal derivation of F(u) = |Vtt| dx, or by letting p — > 1 in the 
definition of the p-Laplacian, p > 1. By analogy with the definition of the first 
eigenvalue of the p-Laplacian on SI, we define the first eigenvalue Ai,n of the 1- 
Laplacian on SI by the minimization problem 



Ai,n = inf / |Vu| dx , 

fue Hi (Si) Jn 

[In \ u \ dx = 1 

where Hl(Sl) is the closure of C^ D (Sl) in the Sobolev space Hl(Sl) of functions in 
L 1 (f2) with one derivative in L . 

The purpose of this paper is the study of the dependence of Ai.n under regular 
perturbations by diffcomorphisms of SI, i.e. we want to compute the first variation, 
the so-called shape derivative, of the functional SI — > Ai,n- General results about the 
stability of A^si under perturbations of SI have been obtained in [17j . In particular 
the authors of [17] found the shape derivative of SI — > Xi.n in the case of regular 
perturbations by diffcomorphisms close to homothctics. We want to extend this 
result to the case of a general perturbation by diffeomorphisms. 

Let us recall some known facts about A^n (see e.g. [ElEQ]). A natural space to 
study A^n is the space BV(Sl) of functions of bounded variations (see, for instance, 
[H [TTJ [T31 [H] ) . By standard properties of BV(Si), we can also define Ai,n by 

Ai,n = inf / |V«|+ / \u\dH n -\ (1) 

\u S BV{Sl) 

1/n \ u \ dx = 1 
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where Vu\ is the total variation of the measure Vu, and iJ™ -1 denotes the {n — 1)- 
dimensional Hausdorf measure. We refer to [8j for a detailed proof of this assertion. 
Note here that if u G BV(SV) and u is the extension of u by in R"\S7, then 
u € BV(R n ) and 

f \Vu\ = f |V«|+ / \u\dH 11 - 1 . (2) 
jk" Jn Jon 

By lower semicontinuity of the total variation and compactness of the embedding 

BV(il) <— > it easily follows from @ that the infimum in ((T|) is attained by 

some nonnegative u G BV(Q). Then u is a solution of the equation Aim = X\^ in 

the sense that there exists a G L°°(f2, R"), ||<7||oo < lj such that 

-div a = X lt Q , 

aWu = |Vu| in fl , and (3) 
{an)u = —u on <9f2 , 

where n is the unit outer normal to d£l, and erVu is the distribution defined by 
integrating by parts J Q (a\7u)v dx when v G Cq 50 (SI) and div a makes sense (eee e.g. 
8, 3J). We then say that u is an eigenfunction for 

We can also express Ai,n in a more geometric way as an isoperimetric type 
problem. We recall that a set C C R™ is said of finite perimeter if its character- 
istic function \c belongs to BV(M. n ). We then define the perimeter \dC\ of C as 
J Rn |Vxc|- Using the coarea formula, we can rewrite {1} as 

Ax,n= _ inf M. (4) 

CCfi, xc£W(I") |0| 

We refer e.g. to [T^ for a proof of this assertion. This infimum is attained by some 
set of finite perimeter, e.g. by a level-set of an extremal for (jT|) , called an eigenset 
or also a Cheeger's set. Note that minimizers for Ai^n touch the boundary dfl 
since, if not, we may blow it up by a factor larger than one, which would decrease 
Ai,r2. Uniqueness and nonuniqueness results of eigensets are in |12[l25j. Concerning 
regularity possible references are [U QUI [HI [TBI ES] • 

To study the variations of A^n with respect to smooth variations of 17, we con- 
sider a smooth vector field V : R™ —* R" and, for small t € R, the diffeomorphisms 
T t defined by 

T t (x)=x + tV(x), (5) 

and eventually the perturbed domains Q t — T t (Q,). We want to compute the 
derivative at t = of the map t — ► Ai^ ■ 

Shape analysis is the subject of an intense research activity. We refer for example 
to Q2] for an introduction to this field. The shape derivative of the first eigenvalue 
A p ,n of the p-Laplacian, p > 1, has been computed in [131 121] : 



^A Pi o t | t=0 = -(p - 1) 



/ 


du p 


Ida. 


dv 



p 

(V, v) dH r ' 



where u p is the unique positive normalized eigenfunctions for X Pt n and v is the unit 
normal vector to dft. What could be the shape derivative of the first eigenvalue of 
the 1-Laplacian is thus not obvious from this formula. 

Our result is the following: 
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Theorem. Let fl be a smooth bounded open subset ofW 1 , V : W 1 — > M" a smooth 
vector-field and fit — Tt(Q), where the Tt's are the diffeomorphisms defined by 0). 
Then 

Ai,n, Ai,n 

as t — » 0. Moreover, if we assume that there exists a unique nonnegative eigen- 
function u € _BV(f2) /or Ai,n such that J Q \u\ dx — 1, i/ien u = |A| xa for some 
eigenset ic 0, and the map t — > X± t n t is differ entiable at t = wii/i 

j /" dTT n ~^ 

j t \u t \t=o = J (divV-(u,DVv) + \i, Q (y,v)) (6) 

where v is the Radon- Nykodym derivative ofVxA with respect to \ (i-£- ^XA — 

^|Vxa| as measures), and d*A denotes the reduced boundary of A, i.e. the part 
of the boundary of A at which can be defined a notion of unit normal vector in a 
measure theoretic sense (see e.g. [2l II H RHI 126] j. 

Let us assume that fi C R™ is strictly convex. We then know from [U [35] that 
there exists a unique Cheeger set A C Cl. Since A has minimum perimeter among 
all the subsets of f2 of finite perimeter and of volume \A\, it follows from [25] that 
dA is C 1,1 . Hence the unit exterior normal vector to dA is defined IF 1-1 a.e. and 
Lipschitz. Its components are thus differentiable at H n ~ l almost every point of 
dA. Note that this vector coincides with — v H n ~ 1 -a.e.. The mean curvature Hqa 
of dA is defined by Hqa = — divg^i/, where Aivqa denotes the tangential derivative 
on dA. We can now write that 

div V - [v, DVv) = div SA V = div g V aA - H aA (V, u), 

where Vqa denotes the tangential part of V, and div ff the divergence operator of 
the manifold (dA, g), g being the metric on dA induced by the Euclidean metric 
(see e.g. [IB]). We can thus rewrite ^ as 

d f dH n ^ 

-7T^i,n t \t=o = / (div g Voa - H dA (V, v) + Ai,n(V, v)) 

at j dA \a\ 

and thus 

d f dTT n ~^ 

-X^^jjX^-HeMu)^-. (7) 

When A = fl, a situation that happens when C R™ is smooth convex and its 
curvature is less that \dQ\/ ((n — l)|f2|) (see [19] when n = 2, and [3] for an arbitrary 
n), formula 10 writes as 

f t Kn t \ t =o = JjXi,n - H dn )(V, v) 

where v is the inner unit normal to dQ. In the particular case when f£ is a ball, 
A = Q and Hqq is constant, so that if we consider a deformation that preserves the 
volume, i.e. a vector-field V such that divy = 0, we get 

d 

Hence a ball is critical for such deformations. 



The following section is devoted to the proof of the theorem. 
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Proof of the theorem 
To simplify the notations, we let A = Ai t n and At = \i.n t . 

Using the change of variable formula for functions of bounded variations [14] , we 
can rewrite At as 



\D(x,t).u v \C(x,t)\Vv\ 

A t = inf ^ j , (8) 

veBV{n) / C(x,t)\v\dx 

where D(x, t) = (DTt(x))^ 1 , C(x,t) = |det (DT t (x))\, and v v is the Radon-Nikodym 
derivative of Vv with respect to |Vu|. Recall that v denotes the extension of v to 
R n by - see ©. As \v v \ = l \Vv\ - a.e, 

\D(x,t)\C(x,t)\Vv\ 

A t < inf J -^— r . (9) 

veBV{n) f C(x,t)\v\dx 



Since \D(x, t)\, C(x, t) — > 1 as t — > uniformly in x G f2, we deduce from © that 

limsupA t <A. (10) 

We let u t G -By (fit) be a nonnegative eigenfunction for At normalized by itt fife = 
1, and v t = u t oT t G BV(Q). Then (w t ) is bounded in BV(W l ). Indeed if we denote 
by T>t the Radon-Nikodym derivative of Vitt with respect to |Vut|, we have 

/ |V«t|= / \(BT^)- x v t \\A^DT^ x \ |VGt|< (l + o(l)) / |V«t| , , 

JR" Jn t Jf2 t (11) 

= (l + o(l))A f , 

and 

v t dx= w t |detDT t _1 |da; = 1 +o(l). (12) 

We can thus assume that the Vt's converge to some nonnegative v G L}{D) in L X {D) 
and a.e. in D, where D is a smooth bounded open subset of K™ containing both 
and the fit's. We denote by u the restriction of v to f2. Then u G SV(f2), u > 0, 
and in view of (fl"2")) . 

it efe = 1 . 



The lower semi-continuity of the total variation and (|11[) then give 

A < J |Vtt| < J |V«| < liminf ^ |V« t | < liminf A t . 
We then deduce with dTUJ) that 

At -> A = / |Vu|, 

as t — » 0, and then that 

/ |V«t| -» A= / |V«|, (13) 
in Jn 
as i — > 0. This proves the first part of the theorem. 
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Let us note for a future use that, as in [17], it follow from (fl3|) that \ Vv t \ — > |Vu| 
weakly in the sense that 



<£|ViJ t |- / 0|Vu| (14) 

for any e C(K n ) with compact support. 

We now prove the differentiability of the map t — > At and the formula ©• We 
use u as a test-function in © to estimate At, so that 



|£>(a;,<)i/|C(a!,<)|Vw| 

A, A ^ ? A, (15) 

C(:r, i)u dec 



where v is the Radon-Nikodym derivative of Vu with respect to | Vm|. Since \v\ = 1 
Vit|-a.e., we can assume that \v\ = 1 everywhere. Direct computations give that 

C{x,t) =det (DT t (x)) = 1 + fdiv7(x) + o(t), (16) 

and 

|D(a;, t)v\ = l-(v; DV{x)v)t + o(t), (17) 
where the o(t) is uniform in x. Thus (|15p becomes 



X + t (divV- (u,DVv)) \Vu\ + o(t) 

At - A < ^ , A 

1+t / udivVdx + o(t) 



Hence 



and 



= t I J divV (\Vu\ - Xudx) -J _(v, DVu)\Vu\J +o(t). 

limsup — < / {{&\vV - {v,DVv))\Vu\- Xu&ivV dx], (18) 

liminf At - A > /{(divV- (v,DVv))\Vu\ - XudivVdx}. (19) 

t^o t Jq 



It remains to prove the opposite inequalities. We use vt as a test-function to 
estimate A, so that 

, , /jV^tl 

X t -X = / |V«t|-A> / \D(x,t)vt\C(x,t)\Vvt\-*2 , 

J ^ Jn / v t dx 

Jn 

where v t denotes the Radon-Nykodym derivative of Vv t with respect to |V«t|. As 
previously, we can assume that \v t \ = 1 everywhere. In view of (|16[) and (jTTJ) , we 
obtain 



|V5t 

A; - A / |Vr,| - / / (divF- (i^.Wwt)) |V«t| - ^ + o(i). (20) 

Vt dx 
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Since divV^ 6 C(O), we get thanks to ([T4|) that 

/ dWV \Vv t \ -> / divV"|Vu| 
Jn Jn 

Independently, 



/ v t dx= I u t \detDT t 1 \dx 
Jn Jn 



with 

det DT^ 1 = 1 - t div V + o(t), 

so that 

/ Vt dx = 1 — t / u t div V + o(t) = 1 —t u div + o(t). 

JO JOt Jo 

Thus 

h - [_\Vv t \+t [ |Vw t | / udivFcfa! + o(t) 



= / |V« t |+*A / udivVdx + o(t), 
Jn Jn 

where the last equality follows from (fT5|) . Hence (j2"0)) becomes 

A 4 -A>i / divF(|Vu| - Xudx) -t / (z/ t , DW t ) |Vut| + o(t). 



Eventually, in view of (|13p and the weak convergence of Vvt to Vu, which follows 
from the L 1 convergence of the Vt to u, we can apply Reshetnyak' theorem [2l [23l [22] 
to get that 

l{x,v t {x))\Vv t \ ->• / g(x,v(x))\Vu\ 
Jn 

for any continuous function g : Ct x S — > R, where S* denotes the unit sphere of R™. 
In particular, 

f(u t ,DVu t )\Wv t \ = [_(v,DVv)\VU\ + o(l). 

JO JO 

Hence 

A* - A > i [ {(divV- (v, DVv))\Vu\ - XudivVdx} + o(t), 
Jn 

and thus 

limsup Xt - A > [{{divV- (v, DVv))\Vu\ - XudivVdx}, 
t^o+ t Jn 

A — A f 

liminf — < / {{divV- (i/, DVv))\Vu\ - XudivVdx}. 

t^o- t J n 

Since by assumption u is the unique normalized eigenfunctions for A, we deduce 
from these two inequalities and (fl8|) . (fT9)) that the map t — > A t is differentiable at 
t = with 

d 



and 



dt 



X t{t=0 = [ {{divV - (v,DVv))\Vu\ - XudivV dx}. 
Jn 
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Eventually, as there always exists an extremal for the problem (|4|), there exists a set 
of finite perimeter A C f2 such that u = |^4| _1 xa- ft then follows from geometric 
measure theory that |Vu| = lA]- 1 !!™'^ (see e.g. [TQ QH US] ) ■ The previous 
formula becomes 

from which we deduce ([6]) using Green' formula for sets of finite perimeter. This 
ends the proof of the theorem. 
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